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Abstract. In this paper we develop an asymptotic analysis for formal and actual solutions of g-difference 
equations, under a regularity assumption, namely the non-collision and non- vanishing of the eigenvalues. In 
particular, evaluations of regular solutions of regular g-difference equations have an exponential growth rate 
which can be computed from the g-difference equation. 

The motivation for the paper comes from a problem in Quantum Topology, the Hyperbolic Volume 
Conjecture, which states that a sequence on Laurent polynomials (the so-called colored Jones function of a 
knot), appropriately evaluated, becomes a sequence of complex numbers that grows exponentially. Moreover, 
the exponential growth rate is proportional to the volume of the knot complement. 

The connection of the Hyperbolic Volume Conjecture with the paper comes from the fact that the colored 
Jones function of a knot is a solution of a g-difference equation, as was proven by TTQ. Le and the author. 
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1. Introduction 

1.1. The goal. The goal of the paper is to intiate an approach to the Hyperbolic Volume Conjecture, via 
asymptotics of solutions of difference equations with a small parameter. The Generalized Volume Conjecture 
links (conjecturally) the (colored) Jones polynomial of a knot to hyperbolic geometry of its complement. 

Since the colored Jones polynomial is a specific solution to a linear g-difference equation, it follows that 
the generalized volume conjecture is the WKB limit of a specific solution of a linear difference equation with 
a small parameter. 

Motivated by this, we study WKB asymptotics of formal and actual solutions of difference equations with 
a small parameter, under certain regularity asymptions. 



1.2. The colored Jones function. A knot in 3-space is a smooth embedding of a circle, considered up to 
isotopy. Two of the simplest knots, the Trefoil (3i) and the Figure Eight (4i) are shown here: 




By the very definition, knots are flexible objects defined up to isotopy, which allows the embedding to 
move in a smooth and arbitrary way as long as it does not cross itself. In algebraic topology, a common 
way of studying knots (and more generally, spaces) is to associate computable numerical invariants (such as 
Euler characteristic, or Homology). Invariants are useful in deciding whether two knots are not the same. It 
is a much harder problem to construct computable invariants that separate knots. 

The invariant that we will consider in this paper is the Jones polynomial of a knot; Q], which is a Laurent 
polynomial with integer coefficients, associated to each knot. The quantum nature of the Jones polynomial 
is apparent both in the original definition of Jones (using Temperley-Lieb algebras) and in the reformulation, 
due to Witten, in terms of the expectation value of a Quantum Field Theory; see PQ IWtj . 

The combinatorics associated to a planar projection of a knot show that the Jones polynomial is a 
computable invariant. However, it is hard to see from this point of view the relation between the Jones 
polynomial and Geometry. In Quantum Field Theory, one often reproduces Geometry by moving carefully 
chosen parameters of the theory to an appropriate limit. 

In our case, we will introduce a new parameter, a natural number which roughly speaking corresponds 
to taking a connected n-parallel of a knot. The resulting invariant is no longer a Laurent polynomial, but 
rather a sequence of Laurent polynomials. 

The colored Jones function of a knot K in 3-space is a sequence of Laurent polynomials 

Jk ■ N — > Zfe*]. 

The first term in the above sequence, Jic(l) is the Jones polynomial of K; see |GL1| . 
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1.3. The Hyperbolic Volume Conjecture. Although knots are flexible objects, Thurston had the idea 
that their complements have a unique decomposition in pieces of unique "crystaline" shape. The shapes in 
question are the 8 different geometries in dimension 3, and the idea in question was termed the "Geometriza- 
tion Conjecture" . The most common of the 8 geometries is Hyperbolic Geometry, that is the existence of a 
complete, finite volume, constant curvature —1 Riemannian metric on knot complements. Thurston proved 
that unless the knot is torus or a satellite, then it carries a unique such metric; see |Th| . 

The Hyperbolic Volume Conjecture (HVC, in short) connects two very different views of knot: namely 
Quantum Field Theory and Riemannian Geometry. The HVC states for every hyperbolic knot K 

lim bgL^nKe^ = 1 yol(5 3_^ 

n— >oo n 2tt 

where vo^S* 3 — K) is the volume of a complete hyperbolic metric in the knot complement S 3 — K . The 
conjecture was formulated in this form by Murakami-Murakami MM following an earlier version due to 
Kashaev, |K|. More generally Gukov (see |Gup formulated a Generalized Hyperbolic Volume Conjecture 
that identifies the limit 

log\J K (n)(e^)\ 

n — >oo ri 

of a hyperbolic knot with known hyperbolic invariants (such as the volume of cone manifolds obtained by 
hyperbolic Dehn filling), for a 6 (0, 1] — Q or a = 1. Actually, the GHVC is stated for complex numbers a. 
For simplicity, we will study asympotics for real a G [0, 1]. 

At present, it is not known whether the limit in the HVC exists, let alone that it can be computed. 
Explicit finite multisum formulas for the colored Jones function of a knot exist; see for example jGLlj . From 
these formulas alone, it is difficult to study the above limit. In a sense, the question is to understand the 
sequence of Laurent polynomials that appears in the HVC. If the sequence is in some sense random, then it 
is hard to expect that the limit exists, or that it can be computed. 

Since the first term of this sequence is the Jones polynomial, and since we know little about the possible 
values of the Jones polynomial, one would expect that there is even less to be said about the colored Jones 
function. 

1.4. g-difference equations. Luckily the colored Jones function behaves in a better way than its first 
term, namely the Jones polynomial. This can be quantified by recent work of TTQ Le and the first author, 
who proved that the colored Jones function of a knot satisfies a g-difference equation. 

In other words, for every knot K there exist rational functions bi(u, v), . . . ,bd(u,v) € Q(u,v) (which of 
course depend on K) such that for all n € N we have: 

d 

Y t W,q)J K {n + j) = Q. 

3=0 

This opens the possibility of studying the g-diffcrence equation rather than one of its solutions, namely 
the colored Jones function. Although the g-difference equation is not unique, it was shown by the first 
author in |Gal| that one can choose a unique g-difference equation, which is a knot invariant. Moreover, 
it was conjectured in |Gal| that the characteristic polynomial of this g-difference equation determines the 
characters of SI_2(C) representations of the knot complement, viewed from the boundary. 

As was explained by the first author on several occasions, asymptotics of solutions of g-difference equations 
would have consequences on the HVC. 

In this introductory article we review the history of asymptotics of solutions of g-difference equations. 

1.5. Asymptotics of differential equations with a parameter. Excellent references for differential 
equations with a parameter are Olver's and Wasau's books; |0] and [Wa . In 1837, Liouville and Green 
independently studied systematically existence of formal (i.e., perturbative) and actual solutions for second 
order differential equations with a parameter; see |Grl IT] . Second order equations are very important for 
classical and quantum physics. 
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In 1908 Birkhoff had the insight to introduce and study arbitrary order differential equation with a 
parameter (see |Blp : 

(1) V {n) + pa n -x(x, p)y (n - 1) + ■■■+ p n a (x, p)y = 

where y = y(x, p), y( n ' means n-th derivative with respect to x (assumed to be restricted to a real interval), 
and p is a large complex parameter, and where the coefficient a,j(x,p) are complex C°° functions with an 
expansion 

oo 

a j( X >p) = ^2 a j,s{x)P~ S 

s=0 

Birkhoff 's working assumption was that the eigenvalues Ai(x), . . . , X n (x) of the characteristic equation 

A" + a n _i :0 (x)A" _1 + • • • + a , (x) = 

were distinct but not necessarily nowhere vanishing. 

In 1926, three theoretical physicists, Wentzel-Krammer-Brillouin studied the second order differential 
equation under the assumption that its eigenvalues do not collide, and developed connection formulas 
linking solutions in the exponential region with those in the oscillatory region. Their method is often referred 
to as the WKB method. 

1.6. Asymptotics of difference equations. As a motivation for our results, let us recall some fundamental 
results of Birkhoff and Trjitzinsky from 1930 on difference equations without a parameter; see jB2j and |BT| . 
A difference equation for a discrete function / : N — » C has the form: 

d 

(2) $>i(n)/(n + j)=0 

3=0 

where aj : N — > C are discrete functions so that ao(n)a.d{n) ^ for all n G N. We will assume the existence 
of asymptotic expansions of a,j(n) around n — > oo for all j = 1, . . . , d: 

a,j(n) ~„_>oo n dj/uJ (a jfi + a^in^ 1/u + a^ 2 n~ 2/uJ + . . . ) 

where u> G N. This certainly holds for ui = 1 if aj are rational functions of n, as is often the case in 
combinatorial problems. 

Due to the nowhere vanishing of ■ ao, it follows that the set of solutions of is a vector space of 
dimension d. 

There are two main problems of difference equations: 

• Existence of formal series solutions tpi , . . . , ipd to . 

• Existence of a basis {tpi, . . . , 4>d} of solutions so that ipk( n ) is asymptotic to V ; fe( n ) f° r large n. 

In |B2j . Birkhoff solved the existence of formal solutions in complete generality (that is, without any 
assumptions on the eigenvalues of the characteristic equation). In a sequel paper |BT| . Birkhoff- Trjitzinsky 
solved the second problem in complete generality. 

Among other things, the formal solutions of Birkhoff lead to the development of differential Galois theory, 
see |vBSj . 

Decades later, the results of Birkhoff and Trjitzinsky on difference equations have found applications to 
cnumerative combinatorics and numerical analysis; see for example Wimp and Zeilberger in |Wil fWZ] and 
references therein. It is not surprising that difference equations are used in numerical analysis, since difference 
equations are numerical schemes of approximating differential equations. In enumerative combinatorics and 
complexity theory, difference equations appear in recursive computation. For example, the number f(n) of 
involutions of {1, 2, . . . , n} (that is, permutations which are a product of 1 and 2-cycles) is given by 

/(n + 2) = /(n + l) + (n + l)/(n) 

with /(l) — 1, /(2) — 2. Using the results of Birkhoff- Trjitzinsky and the fact that /(n) is monotone, it 
follows that 

/w~™^-" /2+ """(' + ^ + ? + ^ + -) 
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for nonzero constants Cj and some K > 0. Actually, the can be computed recursively from the difference 
equation; see [WZ, p. 169]. 

1.7. Asymptotics of difference equations with a parameter. By some historical accident, asymptotics 
of solutions of difference equations with a parameter was not discussed a century ago. The first paper that 
discusses second order difference equation with a parameter appears to be due to Deift-McLaughlin (see 
|DM| ) which was generalized by Costin-Costin to arbitrary order difference equations, jCCj . 

The purpose in this paper is to show that for regular q-difference equations, a regular solution has a 
well-defined and computable exponential growth rate in terms of a relative entropy of the characteristic 
polynomial of the q-difference equation; see Theorem^ below. 

This subject is classical and has been reinvented over the past hundred years by several groups, often 
unaware of each others results. In a sense, the problem of formal solutions of q-difference equations is a 
problem in differential Galois theory; |vPS| . and a problem in numerical analysis; see for example |CCj . 

Our results are hardly new and are contained or can be obtained by minor modifications from results of 
Costin-Costin or from work of Birkhoff and collaborators, |B1I IBTI ICC| . 

Since the presentation in the above papers varies by time and taste, we have decided to give a self-contained 
account of the theory with complete proofs. Hopefully, this will benefit the researchers in Quantum Topology 
and in Analysis. 

1.8. Statement of the results. In this paper, we will describe asymptotics of solutions of q-diffcrence 
equations. 

A q-difference equation for a sequence (/(l), /(2), /(3), . . . ) of smooth functions of q has the form: 

d 

(3) ^^(3 fe ,g)/(fc + i,g) = o 

3=0 

where bj(v,u) are smooth functions and f(k,q) = f(k)(q). 

Before we proceed further, let us remark that q is a variable in @ and not a complex number of absolute 
value less (or more) than 1. In the usual analytic theory of q-difference equations, q is a complex number 
inside or outside the unit disk. 

Moreover, in the GHVC, we need to compute the nth term f(n, q) in the above q-difference equation, and 
then evaluate it at q n — e 27 ™"/™, for a fixed. In other words, in the GHVC, q n is a complex number that 
varies with n in such a way that it stays in the unit circle and approaches 1 as n — > oo. 

With this in mind, e-difference equations (defined below) are obtained from q-difference equations by the 
substitution q = e 2r ™ e where e is a small nonnegative real number, that plays the role of Planck's constant. 

The characteristic polynomial of the q-difference equation @ is 

d 

P(v,\)=Y i b j (v,l)X> 

3=0 

Definition 1.1. We will say that is regular if 

Dsc x P{v, X) ■ b (v,l) ■ b d (v,l) ^ 

for all v e S , where Dsc\P(v, A) is the discriminant of P(v, A), which is a polynomial in the coefficients of 
P(v,X). 

Let Ai(u), . . . , Xd(v) denote the roots of the characteristic polynomial, which we call the eigenvalues of 
(|3*|). It turns out that © is regular iff the eigenvalues X\(v), . . . ,Xd(v) never collide and never vanish, for 
every v G S . Moreover, it follows by the implicit function theorem that the roots are smooth functions of 
v G S 1 . 

Since we are interested in asymptotics of solutions of q-difference equations which, as we shall see, are 
governed by the magnitude of the eigenvalues, we need to partition the circle according to the magnitudes 
of the eigenvalues. 



6 



STAVROS GAROUFALIDIS AND JEFFREY S. GERONIMO 



Let S 1 = Up g -p/p denote a partition of S 1 into a finite union of closed arcs (with nonoverlapping interiors), 
such that the magnitude of the eigenvalues does not change in each arc. In other words, for each p G V, 
there is a permutation o~ p of the set {1, . . . , d} such that 

\K p (i)(v)\ > \K p (2){v)\ > •■■ > \K p (d)(v)\ for all v £ I p . 
The following definition introduces a locally fundamental set of solutions of g-difference equations. 

Definition 1.2. Fix a partition of / as above. A set {ipi, ■ ■ ■ ,ipd} is a locally fundamental set of solutions 
of 121 iff for every solution tp for every p 6 V and for every m = 1, . . . , d there exist smooth functions c v m 
such that 

(4) ^(k,q) = 4(9)ipa p (i){k,q) + --- + (^ d (q)^p(d){k,q) for all (k,q) : q k G I p . 

Theorem 1. Assume that J2J is regular. Then, there exists a locally fundamental set of solutions {ipi, ■ ■ ■ , ipd} 
such that 

• For every m = 1, . . . , d and (k, q) such that q k G I p we have 



n / ka a 
exp -$ m — , - 

a \ n n 



• for some smooth functions <& m with uniform (with respect to x G / = [0, 1]) asymptotic expansion 

oo 
s=0 

where 4> m:S G C°°(T) for all s 

• and leading term 

(5) <p m ,a( x ) = / log(A m (e 27r2 *))dt 

Jo 

where we have chosen a branch for the logarithm of X m . 

Remark 1.3. For every j = 1, . . . , d the smooth functions <fij, s for positive s are uniquely determined from 
the coefficients bj(u,v) of ® by a hierarchy of first-order differential equations along with specified initial 
conditions. On the other hand, the smooth functions $ m are not uniquely determined, since they are obtained 
by a smooth interpolation. Thus, the locally fundamental set of solutions is not uniquely determined from 
the g-differcncc equation, although its asymptotic behavior is. 

It follows from Theorem^ that each locally fundamental solution ijj m (n,q) of the q-difference equation 
© satisfies the GHVC in the sense that for every a G [0, 1] we have: 

lim [\o g \X m (e^)\dt 

n^oo n J Q 

Fix a solution tp of ©■ Theorem ^ expresses ip as a linear combination of f/Vn's in each arc I p . For every 
p G P, let 

(6) Sp = {me{l 4| <4^Q}. 

Later (in Section lfi.3|) we will define the notion of a regular solution to a (/-difference equation. 

As a prototypical example, consider an g-difference equation that satisfies |Ai(i>)| > |A., (i>)| for all j ^ 1 
and all v £ S 1 . Then, any solution that satisfies ci(0) ^ (or more generally, ci has a nonvanishing derivative 
at 0) is regular. 

Remark 1.4. It is possible that S p ^ S p +i. In other words, the restriction of a fixed solution ip to different 
intervals I p may be a linear combination of different ipjS. This is an important phenomenon, referred by the 
name of Stokes phenomenon; see |Waj . 

Our next definition captures the growth rate of regular solutions to regular g-difference equations. 
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Definition 1.5. Fix a collection S = {S p \p e V} of subsets of {1, . . . , d}. The S- entropy 

os ■ [0, 1] — R 

of the g-difference equation © is defined by 

* s (a)= [ \og X s{e 2mat )dt. 
Jo 

where xs ■ [0, 1] — > R is defined by 

Xs(v) = max\\ ap{j) (v)\ if v € I p . 

The entropy of © is the set of functions 

{a s :[0,l]^R | Sc{l,...,d}}. 

Notice that the entropy of a g-difference equation is not a real number, but rather a finite collection of 
functions. 

The main result is the following 

Theorem 2. If f is an S-regular solution of the regular q-dijjerence equation then for every a G [0, 1] 
we have: 

r \og\f(n)(e 2 -^)\ 

iim = as{a.) 

n— >oo 7X 



Finally let us define the J-entropy of a knot. In |Gal| the first author showed that to every knot K one 
can associate a canonical g-difference equation of degree d, and a specific solution of it, namely the colored 
Jones function of K . 

The q-difference equation itself is an invariant of a knot, which (by definition) is determined by the colored 
Jones function of the knot. Thus, any invariant of the g-difference equation is also an invariant of a knot, 
which is determined by the colored Jones function of the knot. 

Definition 1.6. The J-entropy of a knot is the entropy of its associated g-difference equation. We denote 
the J-entropy of a knot K by 

{4* :[0,1]->R | Sc{l,...,d}}. 

1.9. What's next? The paper was written in the spring of 2004. Since then, a number of papers that 
discuss the asymptotics of the colored Jones function have appeared; see |Ga2l IGa31 IGL2I IGL3| . 

1.10. Acknowledgement. The results of this paper were announced in the JAMI 2003 meeting in Johns 
Hopkins. The first author wishes to thank J. Morava for the invitation, and P. Deligne who suggested the 
asymptotic behavior of solutions of g-difference equations. The first author wishes to thank D. Boyd for 
sharing and explaining his unpublished work and also A. Riese, T. Morley, and D. Zeilberger. 

2. e-DIFFERENCE EQUATIONS 

2.1. e-difference equations. In this section, we will translate asymptotics of solutions of g-difference equa- 
tion in terms of asymptotics of solutions of e-difference equations. The latter are defined as follows. 

Fix a positive number e , a compact interval I of R and a natural number d. We will consider functions 
4> : A e0i / — > C with domain 

(7) A toJ := {(fee, e)| k € N, e € (0, e„], ke £ I}. 
Consider the e-difference equation for a function 4> : A £o j — > C 

d 

(8) ^a,(fce,e)0((fc+j)e,e) = 

j=o 

where a, £ C°°(J x [0,e ]). 
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We will assume that for all j = 0, . . . , d, a,j(x, e) has a uniformly (with respect to x) asymptotic expansion 



oc 



(9) a,j(x,e) ~ e ^o ^2aj, s (x)e s 

s=0 

where a hS 6 C°°(I). 

As we mentioned before, e-difference equations are obtained from g-difference equations by the substitution 
q = e 2me where e is a small nonnegative real number, that plays the role of Planck's constant. 
The characteristic polynomial of JHJ is 

d 

P(x, A) = '^2a j {x,0)X J 

3=0 

Definition 2.1. We will say that JSJ regular if 

Dsc\P(x, A) • aa(x, 0) • aa{x, 0)^0 

for all x G I. 

Let Ai(x), . . . , \d{x) denote the roots of the characteristic polynomial, which we call the eigenvalues of 
©• 

It turns out that © is regular iff the eigenvalues Ai(x), . . . , Xd(x) never collide and never vanish, for every 
x G I. Moreover, it follows by the implicit function theorem that the roots are smooth functions of x G I. 

Since we are interested in asymptotics of solutions of e-difference equations which, as we shall see, are 
governed by the magnitude of the eigenvalues, we need to partition the interval I according to the magnitudes 
of the eigenvalues. 

Let / = Upg-p/p denote a partition of I into a finite union of closed intervals (with nonoverlapping 
interiors), such that the magnitude of the eigenvalues does not change in each interval. In other words, for 
each pGf, there is a permutation o~ p of the set {1, . . . , d} such that 

|A CTp (i)(x)| > |A CTp(2) (x)| > ••• > \X ap (d){x)\ for all x E I p . 

The following definition introduces a locally fundamental set of solutions of e-difference equations. 

Definition 2.2. Fix a partition of / as above. A set {ipi, ■ ■ ■ , V^} is a locally fundamental set of solutions 
of © iff for every solution ip : A e j — > C, for every p € V and for every m — 1, . . . , d there exist smooth 
functions cf n e C°° [0, e] such that 

i>(ke, e) = c^(e)V> CTp (i)(fce, e) H ^(e)^ p (d)(^e) for all (fce,e)eA eJ . 

Here, the notation does not indicate the pth power of c m . 

The next theorem summarizes the results of Costin-Costin. 

Theorem 3. ^ |CC| ) Assume that JSJ is regular. Then, there exists a positive e' < £q and a locally funda- 
mental set of solutions {ipi, ■ ■ ■ , ipd} of (|SJ such that 

• For every m — 1, . . . , d and (fee, e) G A e /./ we have 

ip m (ke, e) = exp (e _1 $ m (fce, e)) . 

• for some smooth functions $ m G C°°(I x [0, e']) with uniform (with respect to x G I) asymptotic 
expansion 

oc 

(10) $m(x,e) - e ^o ^</>m, s (a;)e s 

s=0 

where 4> m ,s G C°°(7) for all s 

• and leading term 

(11) exp(</>' (a;)) = A m (x). 
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Fix a solution tb of Theorem [31 expresses ip as a linear combination of the i/>j's in each interval I p . 
For every p £ P, let 

(12) S p = {me{l,...,d}| <4^0}. 

Later (in Section ET2|) we will define the notion of a regular solution to an e-difference equation. 

As a prototypical example, consider an e-difference equation that satisfies |Ai(x)| > |Aj(x)| for all j ^ 1 
and all x £ I — [a, b\. Then, any solution that satisfies ci(a) ^ (or more generally, c\ has a nonvanishing 
derivative at a) is regular. 

Our next definition captures the growth rate of regular solutions to regular e-difference equations. 

Definition 2.3. Fix a collection S = {S p \p G P} of subsets of {1, ... , d}. The S-entropy 

a s : I -> R 

of the e-difference equation JHJ is defined by 

a s{x) = / logxs(t)dt. 
Jo 

where xs '■ I R i s defined by 

Xs(x) =max|A CTj>(j -)(x)| if i£J p . 

Theorem 4. J/ ^ is an S -regular solution to a regular e-difference equation, and x £ I , we have: 

lim elog \ip(x, e)| = a six) 

The next remarks concern the uniqueness of a set of locally fundamental solutions to JHJ • 

Remark 2.4. For every m = 1, . . . , d the smooth functions 4>m,s for positive s are uniquely determined by @ 
and the initial condition (j> m>s (0) = 0. Indeed, applying Taylor series (with respect to e) in JSJ and collecting 
terms, we get for example: 

xl , , l/2C,oW £;=o 0)j 2 A?„(x) + E -=o ^(x, 0)^(1) 



E?=o %(a;,0)jA^(x) 



1 y ^AA P; 

2 Pa PaA 



\\=\ m (x) 



where / x (a:, A) = d/d x f(x, A) and f\(x, A) = d/d\f(x, A). 
Similarly, for s > 1 we have: 

where H s (x) is a function of derivatives of aj{x, 0) and </> mj t for t < s. Notice that the denominator vanishes 
nowhere since the roots do not collide and do not vanish for every x £ I. 

Remark 2.5. If the coefficients a,j(x, e) of the regular e-difference equation © are analytic functions, then 
the functions 4> m .s of Theorem [3] are also analytic, for every m and s. This follows by induction from 
the differential hierarchy which these functions satisfy, and from the fact that the eigenvalues are analytic 
functions. Even though (f) m , s is analytic for every m and s, the series 

00 

^2<p m>s (x)e s 

s=0 

is in general divergent, and the functions $> m . s of Theorem[2]are not analytic. 
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Remark 2.6. Even though the functions <fim,s ar e uniquely determined by the e-difference equation, the 
smooth functions F m (and thus the locally fundamental set of solutions ip m ) are n °t uniquely determined 
by the e-difference equation. The problem is that smooth interpolation is not unique. Recently developed 
ideas of exponentially small corrections might construct a unique set of locally fundamental solutions when 
the coefficients of (JjJJ are analytic functions. We will elaborate on this in a separate occasion. 

2.2. Converting ^-difference equations to e-difference equations. The translation of q-diffcrcncc 
equations to e-difference equations is as follows. If / satisfies the g-difference equation 

d 
3=0 

then set 

q = e 2nl£ , & 3 (e 27r ",e 2T " e ) = aj(x,e) 
and consider the e-difference equation for a function <ft (with domain A eo j for some eo > and / = [0, 2tt]): 

d 

^a i (fce,e)0((fc+i)e,e) = 

3=0 

The following lemma, although elementary, is the key to translating q-difference equations to e-difference 
equations. 

Lemma 2.7. For every (fee, e) E I x [0, eo] we have: 

(13) 0(fce,e) = /(fc,e 2 ™). 

Consequently, for every a £ (0, 1], we have: 



(14) 

Thus, Theorem^ implies The 



lim ilog|/(fc,e 2 " Q/fc )| =a- 1 limelog|0(a,e)|. 



Proof. Observe that aj(ke,e) = bj{e 



27rike c 2ni£ 



), thus (fc,e) — » (f>(ke,e) satisfies the equation 



^6 J (e 2 " fce , e 2 " e )^((fc + j>,e) = 

3=0 

and so does (fc, e) — > f(k, e 2me ). Since solutions with matching initial conditions are unique, i|13|) follows. 
Equation (|14|) follows from equation H13fl by the substitution e — a/k: 

ilog|/(fc,e 27rM / fe )| = ^\og\4>{ke,alk)\=cr 1 e\og\4>{a,e)\. 



□ 



3. Some linear algebra 

In this section we will review some linear algebra. It is obvious that the complex roots of a monic 
polynomial uniquely determine it. It is also known |GLR| that the eigenvalues of a companion matrix 
uniquely determine it, in case they are distinct. 

Consider a companion d by d matrix 



A = 



( 


1 








\ 








1 






























1 




\-ao 


-ai 


-a 2 ■ ■ 


• -dd- 


J 
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The characteristic polynomial of A is 

d-i 

A d +$>^ 

3=0 

with roots Ai, . . . , A^- Let M = (A* _1 )i J - be the Vandermonde matrix, and D = diag(Ai, . . . , Xd) be the 
diagonal matrix with diagonal entries Ai, . . . , Xd- 

Lemma 3.1. If a companion matrix has distinct eigenvalues, then with the above notation we have: 

A = MDM' 1 

Proof. Observe that Vj — (1, Xj, . . . , X d -^ 1 ) T is an eigenvector of A with eigenvalue Xj. Thus, M = (vi , . . . , V4) 
and AM = MD. The result follows. □ 

Now, consider a companion matrix A(u) whose entries in the bottom row are smooth functions in a 
variable u, with roots Ai(m), . . . , Xd(u) which we assume are distinct for all u. 

The next lemma is a key estimate for the norm of long products of slowly varying matrices. In the 
language of physics, A(u) is the transfer matrix and A(n) . . . A(2)A(X) is the transition matrix. 

Lemma 3.2. Assume that the eigenvalues Ai (it), . . . , Xd(u) of A{u) are distinct for all u and 

max sup |Aj(u)| < 1 + Ce. 

3 u 

Then for m < n, ne € I , we have 

\\A(n)A(n - 1) . . . A(m)\\ < C 

Proof. By Lemma II we have: 

A(u) = M{u)D{u)M{u)~ l 

If m < n, it follows that 

A(n)A(n - 1) . . . A(m) = M(n)D(n)M(n)~ 1 M(n - l)D(n - l)M(n - l)" 1 . . . M{m)D{m)M(m)~ 1 , 

which implies that 

\\A{ri)A{n-l)...A(m)\\ < \\M(n)\\M(m)~ l \\ ■ \\D{n)\\ . . . \\D{m)\\ 

||M(n) _1 M(n - 1)|| . . . \\M(m + l) _1 M(m)||. 

Now, we have 

\\D(k)\\ < 1 + Ce for k = m,...,n 
\\M(k)M(k- ly 1 ]] < 1 + C'e by Lemma 1331 

If / = [a, b], using the fact that ne, me £ /, we obtain: 

\\A(n)A(n - 1) . . . A(m)\\ < (1 + Ce)"-" l (l + C'e)"-" 1 

< (1 + C"e) 2[n ~ m) 

< (1 + C"e) 2 ^ 

< e 2C"(6-a)_ 

□ 

Lemma 6.6. IfM= (x l f 1 ) lJ and N = (y) are Vandermonde matrices, such that M is nonsingular, 

then 

V3 ~ x i 
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4. Existence of formal solutions 
In this section we will prove that JHJ has a unique set of formal solutions. Let us define those first. 
Definition 4.1. A formal series ip(x,e) is one of the form 

(15) Vi(x,e) = exp(e- 1 £^( a ;)e s ) 



s=0 



where <j) s £ C°°(I) are smooth functions for all s. 



Note that e\ogip(x,e) lies in the ring C°°(/)[[e]] of formal power series with coefficients smooth functions 
on /. 

Note further that if ip(x, e) is a formal series, so is ip(x + je, e) for every j £ Z, where the latter may be 
defined using the Taylor series of (f> s (x + je) = H^* O^)-?'* 6 *- 

= f{x, e) exp ^{x)j + ^( X )j + e + + + ^Mj 3 ) £ 2 + . 

It follows that if ^(x, e) is a formal series, then the ratio ip(x + e, e)/ip(x, e) lies in the ring C°°(/)[[e]]. 
Using the language of difference Galois theory (see |vPSl p. 4]) this implies that 

Lemma 4.2. C°°(/)[[e]] is a finite difference ring, under the map x — » x + e. 

Definition 4.3. We say that a formal series ip of l|15fl is a formal series solution to (0 iff 

d 

(16) 77 r E °i(*> e * + J e > £ ) = e ^WM}- 

It is easy to see that if tp is a formal solution to JSJ, then the leading term cf>o satisfies the equation 

(17) exp^x)) = A(z) 
where X(x) is an eigenvalue of JHJ). 

Proposition 4.4. Assume that (jSJ) is regular. Then, (|SJ) /ias d unique formal series solutions tpi, . . . ,ipd 

with leading terms corresponding to the eigenvalues of JSJ). 

Proof. First we need to show that l|l(j|) is indeed an equation in the power series ring C°°(/)[[e]], i.e., that 
the terms involving negative powers of e cancel. 

Suppose that ip is given by (|15fl . It follows from the calculation preceding Lemma l4~21 that for every s £ IN, 
we have: 

(18) coes(e^^ +J£ > e) )={ eMmx)) ^ S = ° 

\ ip(x,e) J [jexp(j(j)' (x))(j)' s (x) + terms s if s>0 

where coeff(e s , <?(e)) denotes the coefficient of e s in a power series g(e), and where terms s is a polynomial in 
the derivatives of <f>t for t < s. 

Expand the terms of Equation (jl6|) into power series in e using the above equation and 10 , and collect 
terms of powers of e. It follows that Hlfi[) is equivalent to a hierarchy of first order differential equations: 

d 

^2aj(x,0)exp(j(f>' (x)) = 
3=0 

d 

aj(x, 0)j exp(j(j>' (x))(f)' s (x) + Terms 8 = 

3=0 
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where for positive s, Terms s is a polynomial in the derivatives of 4>t and aj(x,0) for t < s. 

Now fix an m G {1, . . . , d}, and choose 4>m,o such that exp((f>' m (x)) = X m (x), where X\(x), . . . , \d(x) are 
the eigenvalues of (jHJ- Since (JHJ) is regular, it follows that the roots X\(x), . . . , Xd(x) of the characteristic 
polynomial P{x, X) never collide, and never vanish for x G /. Thus, 

? (\±P(x,\)) = J2j ai (x,0)\l(x) 

\ aA /A=A m (x) j=0 

for all x E I. Thus, after we choose 4> m ,0i it follows that we can find functions m , s for s > that satisfy 
the above hierarchy. Moreover, for every m, the sequence of functions 4> m , s is uniquely determined by the 
above hierarchy and the initial conditions 4> m .s(0) = 0. □ 

4.1. An alternative formal series. In this section we present an alternative, and slightly more general 
form, of formal series. In case of regular e-diffcrcncc equations this alternative form will not be needed. 
However, when eigenvalues collide or vanish, one must use this alternative form of formal series. Thus, in 
the present paper we will not use this alternative form of formal series, and the reader may skip this section. 

Definition 4.5. An alternative formal series ip(x, e ) 1S one of the form 

oc 

(19) ^(x, e )=exp(e~ 1 0(x))^^(a ; ) e s 

s=0 

where (f>, <f> s G C°°(I) are smooth functions for all s, and 4>o{x) ^ for all x E I. 

In the remainder of this subsection, we will refer to alternative formal series simply by formal series. 
Note that if ip(x, e) is a formal series, so is ip(x + je, e) for any j G Z, where the latter may be defined 
using the Taylor series of 4> s {x + je) and (f)(x + je) around x. It follows that 

$(x + je, e) = exp (e -1 0(a:)) (0o(a) + (<t> {1) {x)M x ) + ffi ( x ) + M x )) e + ■•■)■ 

Moreover, if tjj(x, e) is a formal series, then the ratio ip(x + e, e)/ip(x, e) lies in the ring C°°(J)[[e]]. This 
follows from i|2U|) and the following computation, valid for every j G N: 

4>(x+je,e) = exp{e- 1 <t){x+je))(Mx+je) + <l)i{x+je)e + 0{e 2 )) 



exp 



(e-^ix) + ^{x)j + ^-^fe + 0(e 2 )) (c/> Q ( x ) + ( x )j + M^))e + 0(e 2 ) 



= Hx, e) exp(^) (x)j) (l + [^j 2 + ^ j e + 0(e 2 )j 

In analogy with Lemma 14.21 this implies that 

Lemma 4.6. C°°(/)[[e]] is a finite difference ring, under the map x — » x + e. 

Definition 4.7. We say that a formal series ip of Q19[l is a formal series solution to (JHJ) iff 

d 

(20) -j- ]T aj (x, e)${x + je, e) = G C°°(I)[[e]]. 

It is easy to see that if tp is a formal solution to JSJ), then the leading term <fi satisfies the equation 

(21) exp(<f>'(x)) = \{x) 

where X(x) is an eigenvalue of JHJ- 

In analogy with Proposition 14.41 we have the following: 

Proposition 4.8. Assume that (jHJ) is regular. Then, (JHJ) has d unique formal series solutions tpi, . . . ,ipd 

with leading terms corresponding to the eigenvalues of J^J. 
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5. Proof of Theorem El 

In this section we prove Theorem [21 The strategy is to 

(a) prove that there exists a solution tpi with the stated properties where \i{x) is an eigenvalue with 
maximum magnitude. 

(b) use this solution tpi to reduce Theorem to the case of a e-difference equation of degree one less 
than the original one. 

(c) prove that the constructed set of solutions is a locally fundamental set. 

Without loss of generality, we will assume that the eigenvalues of (jHl satisfy the inequality: 

|Ai(x)|>|A 2 (x)|>-.->|A d (x)| 
for all x E I. Otherwise, we can partition / into subintervals where this is true. 

5.1. Existence of a solution corresponding to the eigenvalue of maximum magnitude. Consider 
first a formal solution tpi of © given in Proposition ^. 41 which satisfies (|15|l and (|1 T|> . Consider the smooth 
functions G C°°(I) of fT7)l. 

The proof of the following lemma (due to Borel in case </>i iS are constant functions, for all s) can be found 
in Lemma 2.5]: 

Lemma 5.1. There exists a smooth function $i G C°°(I x [0, eo]) such that we have (uniformly in x G I): 

oo 

$l(x,e) - e ^ ^2<j) liS (x)e s . 
j=o 

Now, consider the unique solution ipi of JSJ with initial conditions 

^i(fce, e) = exp(e _1 $i(fce, e)) for k = 0, . . . , d - 1 

and for small enough e, where without loss of generality, we assume that / = [0, b]. 

Of course, for large k it may not be true that tpi(ke, e) = exp(e _1( l > i(fce, e)). The next proposition estimates 
the error, uniformly with respect to k: 

Proposition 5.2. There exists an e' > and constants C s such that for all (fee, e) £ A e <j and all s G INI, 
we have (uniformly in k): 

V>i(fce,e) 



1 



< C,e s 



(22) 

exp(e _1< l'i(fce, e)) 
Proof. Let us make a change of variables: 

(23) 6=* 

ipi 

where 

4>i(x,e) = exp ^e _1 $i(a;,e)^ . 

We will show that for a fixed sq, and for every s there exists a constant C s such that for all (fee, e) G A eo j 
we have: 

(24) \9(ke,e)-l\<C s e s+1 - S0 , 
Since satisfies it follows that 9 satisfies 

d 

(25) ^(fce,e)0((fc+j)e,e) = O 
where 

t / \ / x"0i(^ + J'e,e) 
6j(x,e) =a j (a;,e) — f - — . 

ipi{x,e) 

It is easy to see that 
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• bj(x,e) <E C°°(I x [0, Co]), has uniform (with respect to x) e-asymptotic expansion as in 0, 

• b jiS (x,0) = aj(x,0)X{(x), 

• and since ip is a formal series solution to JSJ) and <3?i is given by Lemma |5~T1 if follows that for every 
s we have: 

d 

(26) 5>(a;,e)=0(0 

i=o 

The characteristic polynomial of (|25[1 has roots fj, m (x) := A m (a:)/Ai(a;) for j = 1, 
so is <|25ll . Notice that Ai(a;) vanishes nowhere since (|25H is regular. 

We now show (|24() . Let us write the difference equation l|25|l in matrix form 

(27) 

where 



, d. If JHl is regular, 



Q((k+ l)e,e) = A(fce,e)6(fce,e) 



6»(a; + e,e) 
6(x + 2e,e) 



A(x,e) 



( o 






\9(x + (d-l)e,e)J 



-c d -i(x,e)J 



\~c (x,e) -c%(x, e) -c 2 (a;,e) 

and Cj(a;, e) = frj(x, e)/bd(x, e). Iterating, we obtain that 

6(fce, e) = A(ke, e)A((k - l)e, e) . . . A(e, e)9(0, e) 

for all k > 1, where 8(0, e) = 1, a column vector with all entries equal to 1. 
Equation (|26() gives: 

(28) A(fce,e)l = l + e s E fe (e) 

where ||Efc(e)|| < C uniformly in k and e. Feeding in the above equation, we obtain: 

fc-i 

(29) Q(ke, e) = 1 + e s E fc (e) + A(fce, e)A((fc - l)e, e) . . . + l)e, e)E J (e). 

i=i 

Now, let us look at the roots fi±(x,e), . . . , Hd{%, e) of 

d 

i=o 

Since max.,- sup^gj |^-(:r,0)| = 1, it follows that 

max sup \[ij(x, e)| < 1 + Ce. 

Since ke lies in /, a compact set, Lemma fc>. 21 and Equation (|29|l imply that 

(30) ||O fc (e) - 1|| < kC' s e s+1 < C s e s - Sa+1 

for all ke £ /, where So — 1. This completes the proof of l|24|l . 

Equations l|24|l and 11' I'D differ in the presence of sq . It is easy to see that if / is a function such that for 
a fixed sq and any s £ N we have: 

s+s 



then 



|/( £ )-^C t£ t |<^+ 1 , 



|/(e)-^C t e 4 | < (D s + \C s+1 \ + ...\C s+So \y 



This observation shows that i|24|) implies (|22|l and concludes the proof of the Proposition. 



□ 
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Proposition 5.3. (a) There exists a smooth function $1 G C°°(J x [0, e']) such that 

(a) For all (ke,e) G A e </, we have: 

ipi(ke, e) = exp (e^^fee, e)) , 

(b) $i has an asymptotic expansion (uniform with respect to x): 

$1(2;, e) ~ e ^ exp ^e" 1 ^ M (x)e s ^ 

where (pi tS are as in Lemma \5.1\ As a result, we have an asymptotic expansion (uniform with respect 
to k): 

^l(ke, e) ~ e ^ exp ^e _1 ^ M (fce)e s ^ . 

Proof. Consider the change of variables 9 as in l|23|) ■ 

Due to our choice of initial conditions it follows that for every fixed k = 0, ...,d — 1, the function 
e — > 9(ke, e) is smooth. Using this and the smoothness of the coefficients of JSJl, it follows that for every fixed 
fc, the function e — > 9(ke, e) is smooth. 

So far, the function 9 is defined on A^f S T : 



(0,e) 



fc 

x=0 


x=e 




x=2e 




x=3e 




x=4e 




X 



which is a union of line segments in a rectangle / x [0, e'], and satisfies l|22[) . 

The complement of these line segments in [0, e'] consists of an infinite union of open triangles, together 
with the horizontal segment 7x0. We can smoothly interpolate 9 inside these open triangles so that it is 
defined on I x (0, e'] and 

(31) \9(x,e)-l\<C s e s 

for all (x, e) G I x (0, e'] and all seN. 

Let us extend 9 to I x [0, e'] by defining 9(x, 0) = 1 for all x G I. 

We claim that 9 is smooth on I x [0, e']. We need only to check this at the points (x, 0) for x G I. 
This follows easily from (|3i|) . For example, to check continuity at (x,0), consider a sequence {x n ,e n ) such 
that lim, woo (x„, e„) = (x,0). Then, (|3i|l for s = 1 implies that \9(x n ,e n ) — 1| < C s e„ ~ n -»oo 0, thus is 
continuous at (x, 0). Using (|31(l for s + 1 it follows that d s /de s 9(x, e)| e =o = for all s > 0, and we find that 
9 has an e-asymptotic expansion (uniform with respect to x): 

6{x, e) ~ e _,.o 1 

Restricting further e' if needed, we may assume that \0(x, e)| > for all (x, e) G I x [0, e']; in other words 
log#(x, e) makes sense for all (x, e) G I x [0, e']. 

Now, we can finish the proof of Proposition |^| 

Let us define 

$i (x, e) = <t>i (:r, e) + e log0(x, e) 
Then, 123J1 implies (a). Since 6(x, e) is asymptotic to 1 (uniformly on x), it follows that &i(x, e) is asymptotic 
to $i(ir,e). Using the asymptotic of <SE>i given by Lemma r5.il (b) follows. □ 
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5.2. A reduction to an e-difference equation of smaller degree. We will now prove Theorem [3] by 
induction on the degree d of the e-difference equation. For d — 1, it follows from Proposition 15.31 The 
inductive step is the next Proposition. 

Proposition 5.4. Assume that Theorem holds for regular e-difference equations of degree less than d. 
Then it holds for regular e-difference equations of degree d. 

Proof. Consider a e-difference equation (JHJ) of degree d. We will use the solution tpi of it constructed in 
Proposition 15.31 to reduce it to an equivalent equation of degree d — 1, and an inhomogeneous e-difference 
equation of degree 1. 

Consider the dependent change of variables 

(32) e = -f 

This is well-defined since ipx is nowhere zero. Then, cj> satisfies JHJ iff satisfies 

d 

(33) ^6 i (fe e) e)^((fe+i)e,e)=0 



where 



bj(x, e) = dj(x, e)- 



ipi(x,e) 

The characteristic polynomials of (JSJ) and l|33|l are related by 

= Ai(x) d P il (A/A 1 ( a ;)). 

As in the proof of Proposition I5.3I it is easy to see that 1)33(1 is a regular e-difference equation. Moreover, it 
is easy to see that Theorem |3 holds for JHJ iff it holds for 1)33(1 . Indeed, check that the change of variables 
given by l(32|l preserves the asymptotics of the solutions of |JSJ and ((33|) . 

Thus, it suffices to work with ((33|) . In that case, 8 = 1 is a solution of l(33|) . since tp is a solution of JSJ- 
It follows that 

d 

(34) ^(z,e) = 0. 

i=o 

(Compare this with 1(26(1 ~). Let us define 

(35) C(fcc, e) = 0((fc + l)e, e) - 9(ke, e). 
Then, we get that C is a solution of the e-diffcrence equation 

d-i 

(36) ^ Cj (fce,e)C((fc+J>,e)=0 

3=0 

where 

d 

c s (x,e) = bj(x,e). 

j=s+l 

The characteristic equations of ((33(1 and l(36|) are related by 

d d~l 

^6 j (x,0)A^ = (A-l)^c J (x,0)A^ 

Since co(x, e) = X)j=i e ) = — °o{x, e) (by H34|) 1 and Cd(x, e) = 6d(^, e), the same arguments of Proposition 
IO imply that igU) is regular, assuming that is regular. 

By the induction hypothesis, it follows that ((36(1 satisfies Theorem [21 
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For the remainder of this section, fix a solution ( of (|36[) which satisfies the properties of Theorem [31 
In other words, £ satisfies (|33|) and C,(ke,e) = exp(e _1 Z(fce, e)) where Z is a smooth function with uniform 
(with respect to x) asymptotics: 

oo 

Z(x, e) ~ e ^ ^2 Z s (x)e s . 

Lemma 5.5. There exists a formal solution 

9{x,e) = exp(e~ 1 9(a;,e)) 

of <|36[1 such that 

oo 

&(x,e)=Y / ®s(x)e s . 

s=0 

Proof. We need to solve the formal power series Equation 

(oo \ / oo \ / oo \ 

e- 1 Q s(x + e)e s - exp e" 1 ^ G s {x)e s = exp e" 1 Z s {x)e s 
s=0 / \ s=0 / \ s=0 / 

for O in terms of Z. Using the Taylor expansion <do(x + e) = Qq(x) + <d' (x)e + 0(e 2 ) it is easy to see that 
the above equation equals to 

exp (e^eoO) + 0(1)) - exp (e^Go^) + 0(1)) = exp (e- 1 Z {x) + 0(1)) 

from which follows that 0o = Zq. Dividing the equation by exp(e _1 0o) we get an equation in formal power 
series with nonnegative powers of e. Moreover, the coefficient of e s in that power series (for s > 0) equals to 

(exp(ej,(i)) - l)exp(e s+1 (x))H s (x) 

where H s (x) is a function of Qj and Zj for j = 1, . . . , s. Since exp(9g(x)) = exp(Z' (x)) is an eigenvalue of 
it is never equal to 1. 

This and induction prove the lemma. □ 

Lemma 5.6. (a) There exists a solution to the equation 

(37) ((ke,e) = 9((k + l)e,e)-d{ke,e) 

for 9 in terms of £ with appropriate initial condition, 
(b) For all (fee, e) £ A^, 7 we have: 

£(ke,e) ~ c _>o exp(Z(fce, e)) 

where 

exp(Z(a;. e)) = exp(0(a; + e, e)) — exp(0(x, e)) 

Proof. Fix e > and let / = [a, b}. Consider a natural number k such that ke £ I and (k + l)e £ I. This is 
equivalent to k\ < k < k 2 where k\ and k 2 are natural numbers that depend on e and /, although we do not 
indicate this in our notation. 

Then equation l|37[l implies that 

C(fc l£ ,e) = 9{{k 1 + l)e,e)-9{k 1 e,e) 

C((fci + l)e,e) = 6((k 1 + 2)e,e)-6((ki + l)e,e) 

C((fc 2 -l)e,e) = 9(k 2 e,e)-9((k 2 -l)e,e) 

Summing up, we obtain that 

fc 2 -i 

8(k 2 e,e) = 9(k 1 e,e)+J2 C(je,e). 

Choose initial conditions so that 9(kit,e) — exp(Z(fcie, e)). This completes part (a). 

Part (b) follows by a telescoping calculation. □ 
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To finish the proof of Proposition l5.4l it suffices to show that the solution £ of Lemma \b. 61 is asymptotic 
to the formal solution of Lemma 1531 
Since 

C(ke, e) ~ c ^o exp(Z(fce, e)) 

and 

exp(Z(fce, e)) = exp(9((fc + l)e, e)) - exp(9(fce, e)) 
it follows by the definition of 9 given in Lemma 15 . 61 and by a telescopic sum, that: 



fc-i 

6{ke,e) = fl(Aie,e)+X)C(7'e,e) 

k-l 

~ e ^ Oihe, e) + J2 (exp(6((j + l)e, e)) - exp(9(je, e))) 

= 6»(fcie, e) + exp(9(fce, e)) - exp(9(£e, e)) 
= cxp(9(fce, e)) 

This concludes the proof of Proposition EH □ 

5.3. The solutions form a locally fundamental set. Let us summarize what we have obtained so far. 

Consider a partition I = U pe pl p of I = [xq, xp] given by I p — [x p , x p +i] for p — 1, . . . , P — 1, and consider 
a permutation a v of {1, . . . , d} such that 

\K p (i)(x)\ > \K p (2)(x)\ > • • • > \X ap(d) {x)\ for all x e I p . 

We have constructed solutions smooth functions $ m for m = 1, . . . , d with asymptotic expansion given by 
(|TU|) and {nj. 
Let us define 

(38) ip m (x,e) = exp (e _1 $ m (a;, e)) , 

where $ m are smooth functions with asymptotic expansions as in Equations l|10[) and 

Moreover, we have shown that for every interval I p (as in the discussion prior to Theorem 0), 

{ipi(ke,e), . . .,ipa(ke,e)} 

is a set of solutions of JHJ) when ke E I p . 

Fix a solution ip(ke, e) of (JHJ) and an interval / p . The following lemma certainly implies that {ipi, ■ ■ . , ipd} 
is a locally fundamental set of solutions of JSJl. This concludes the proof of Theorem [3J □ 

In addition, the next lemma motivates the definition of a regular solution, given in the following section. 

Lemma 5.7. (a) Fix ip and I p as above. Then, there exist smooth functions such that 

(39) ip(ke,e) = (e)Vv p (i) (ke, e) + • ■ ■ + c p d (e)^ lTp{d) (ke,e) 
for all ke £ I p . 

(b) Moreover, for every p and m = 1, . . . , d we have 

(401 c p (e) = ^p-^ m ^ x P e ~ e ' e \ P ( e ) 

for some smooth functions 7^, with the understanding that Vv_i(m) — 1- Here [x] is the largest integer 
smaller than x, and 

xi= [^1+1. 
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Proof. Without loss of generality, let us assume that a p (j) — j for j = 1, . . . , d. Equation i|39|) is a linear 
equation in c^, with solutions 

det W m (x e p ,e) 



where I p = [x p ,x p+1 ], 



W{x,e) 



and 



ipi(x + e,e) 



W m (x,e) 



( ipi(x,e) 
ipl (x + e, e) 



detW(x e ,e) 



i> m (x,e) 
ip m (x + e,e) 



\Vi(ar+(d-l)e,e) ... ip m (x + (d - l)e, e) 



V>(ar,e) 
V»(a; + e,e) 

V»(a; + (d - l)e, e) 



\tpi(x+(d-l)e,e) . 

where the ip's are in the roth column of W m . 

We will show that for small enough e, W(x, e) is nonsingular. 
Using Equations (JSHl, 03 and JTTJ, it follows that 



i>d{x,e) \ 
V»d(a; + e,e) 

ifi d (x+(d-l)e,e)J 

V>d(x + e, e) 
ipd(x+(d-i)e,e)J 



— r — -det(W^)) = det(exp(j^(x))) + 0(e) 

^i(ar,e) . . .ipd(x,e) 

= det(X m { X y) + 0(e) 

= ±J[(X i (x)-X j (x))+0(e) 

uniformly in a;, where A m are the eigenvalues of (jHJ. Since (JSJ is regular, its eigenvalues never collide. 
Similarly, using Equation l|38|) . we have: 



^i(x^e, e) . . . ^ m (a£e, e) . . . Vd(^e, e) 



det(W m (x;e, e)) = det B m . (x e p e, e) + 0(e) 



where 



Ai(x) 1 ^ 



^(x,e) 
^i(a; + e,e) 



Ad(x) 1 ^ 



\Ai(a 



^(z + (d - l)e, e) ... \d{x) d -^) 



where the ^'s are in the roth column of B m j. Thus, 

1 det(B mfi (x p e, e)) 



(41) 



^ m (a£e,e) n j#m Aj(^e) - X m (x p e) 



0(e). 



The idea now is to move the recursion relation backwards d times. Using the solution ij)(ke, e) for ke G I p -i 
will allow us to compute the smooth functions c^. 

In detail, consider the matrix B m ^(x p e, e) and move the recursion relation backwards once. Using Equa- 
tion l|27|l and the fact that the the jth column of B m j(x p , e) for j ^ ro is an eigenvector of A(x p e, e) (up to 
0(e)) with eigenvalue Xj(x p e), it follows that 

B mfi (x p e, e) = A(x e p e, e)B mA (x e p e - e, e) + 0(e). 

Iterating d — 1 more times, it follows that 

B m ,o(x £ p e, e) = A(x f p e, e) . . . A(x £ p e - (d - l)e, e)B mA (x p e - de, e) + 0(e). 
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Since x e p e — de G I p -i, it follows (as in the computation of W m (x,e) above) that: 



det(B m (x e p e - de,e)) 



det 



V' CTp _ 1 ( m )(a;pe-e,e) 
(\ 1 {x p e-e) 1 - d . 
\ 1 ( Xp e-e) 2 - d . 



y%-i(^)K £ - de ' £ ) 

r p _i(m)Ke+(l-rf)e,e) 



i-d 



\Ai(a£e-e 
J] (Ari^e-^-AT 1 ^-^) 

This, together with Equation (|4*T)l proves (|i*U|> 



A^e-e) 1 "^ 
A d (x;e-e) 2 - d 

A d (a£e - e) d - d ) 



□ 



6. Regular solutions and their asymptotics 
In this section we discuss regular solutions of q and e-difference equations and their asymptotics. 

6.1. Regular solutions to e-difference equations. According to Lemma f5. 71 a solution ip to (JSJ deter- 
mines a collection 5 = {S p \p G V} of subsets of {1, ... where 

S p = {me{l,...,d}| T :/0}. 

Definition 6.1. Fix a collection S = {S p \p £ P} of subsets of {1, . . . ,d}. We say that a solution ip of a 
regular equation ijSJ is S-regular iff for every p£ P we have: 

. 7 p = if m g S p . 

• 7^ are not flat at for all m £ S p . That is, some derivative of 7^ at e = does not vanish. 

• For every p £ V there exists an element rj (p) G S p such that 

W(p)( x )\ > \^j( x )\ f° r au " 3 ^ Sp — {v(p)}i x S Interior(Jp). 
In other words, in the interior of the interval I p , and among the eigenvalues Xj(x) for j G S p , there 
is a unique eigenvalue of strictly maximum magnitude. 
We will say that a solution to © is regular if it is ^-regular for some S. 

6.2. Asymptotics of regular solutions of e-difference equations. 

Proof, (of Theorem 0J) Let ip be an ^-regular solution to JSJ. Let us assume that S — {1, ...,d}, and 
|Ai(x)| > |A m (x)| for m = 2, . . . , d and all x in the interior Interior(I) of the closed interval / = [0, b]. Fix 
an x £ I. 

Then, we have: 

ip(x, e) = ci(e)ipx(x, e) H h c d {e)ip d {x, e). 

for a; = fee, where ci(e) = cie™ 1 + 0(e™ 1+1 ), and ci 7^ 0. 
Then, we have: 



(42) 

Recall from Theorem [21 that 



ip{x,e) = c 1 (e)ip 1 (x,e) 1 + 



\ - Cj(e) ip m {x,e) 



1p m (x,e) 



exp(e _1 $ m (a;,e)) 
log A m (a:). 
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Thus, 

Rfi($m,o)'(aO = Re(logA m (x)) = log|A m (x)|. 

Combined with |A m (x)| < Ai(x) for in > 2, and < J > m ,o(0) = = < I>i i o(0), it follows that 

Re($ m>0 )(x) < Re($i, )(x) 

for all x € I and ^ 

Re($ 1 , )(x) = / log|Ai(t)|dt. 
Jo 



Therefore, 



and 



lim fl + V c J"( e ) ^m(x,e) \ = ^ 



lim elog |ci(e)| = 0. 



Thus, Equation (|42[l implies that 



lim elog \ip(x, e)| — lim elog \ipi(x, e)| 



lim elog|exp(e $ m (x, e))| 

e— +0+ 

lim elog|exp(e- 1 $ m (x,0) + O(l))| 
Re($i, (x)) 
Iog|Ai(t)|dt 
The result follows. 

In the general case, we partition the interval I p as in the discussion prior to Theorem [3] and repeat the 
above proof using Equation (|43[) . The result follows. □ 

6.3. Asymptotics of regular solutions of g-difference equations. First, we need to define what is a 
regular solution to a g-difference equation. 

Consider a solution ip of a ^-difference equation and a partition of S 1 as in Section 11.81 Then, at each 
interval I p , we can write the solution as a linear combination of fundamental solutions, as in Equation Q. 
Let S p be the indexing set of the fundamental solutions that we use in each interval I p ; see (O- 

Suppose that the partition of S 1 is given by I p — [e 2nlXp , e 27rlx p+ 1 ] for p — 0, . . . , P — 1, with qo = 1, qp = 

Then, with q = e 2rr2e it turns out that for every p and m = 1, . . . , d we have 

(43) c* («) = \ {m){ fr £ : t] im 

for some smooth functions 7^, with the understanding that Vv_i(m) = 1- 

Definition 6.2. Fix a collection S — {S p \p £ V} of subsets of {1, . . . ,d}. We say that a solution ip of a 
regular equation J3Jl is S-regular iff for every p £ V we have: 
. 7*=0ifm£S p . 

• 7^ are raoi /?ai at for all m £ S p . That is, some derivative of 7^ at e = does not vanish. 

• For every p £ V there exists an element rj(p) £ S p such that 

|A 7 ( p )(u)| > |Aj(u)| for all j £ S p - {7(f) }, v £ Interior(J p ). 

In other words, in the interior of the interval I p , and among the eigenvalues Xj(v) for j £ S p , there 
is a unique eigenvalue of strictly maximum magnitude. 
We will say that a solution to © is regular if it is S'-regular for some S. 

Proof, (of Theorem |2l It follows from Equation i|14|) of Lemma 12. 71 and Theorem 0] □ 
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7. Applications to Quantum Topology 

7.1. The A-polynomial of a knot and its noncommutative version. In this section we discuss general 
features of g-difference equations for the colored Jones function of a knot. 

The coefficients of the q- difference equations are rational functions of q and q n = Q. In order to simplify 
the typesetting, we will give the q-difference equation 

d 

5>(g n ,g)/(n + j) = 
j=o 

in operator form 

(44) Pf = 
where 

d 

P = Y i b i {Q,q)E? 
i=o 

and where the operators E, Q and q, act on a discrete function / : N — > Z[g ± ] by 

(?/)(n) = qf(n) (Qf)(n) = q n f{n) (Ef)(n) = f(n + 1). 

Note that q commutes with Q and E, and that EQ = qQE. 

It follows by definition that the characteristic polynomial chP(t>, A) of (|44|) is obtained from P by setting 
9=1, replacing (E,Q) by (X,v). In other words, we have: 

d 

chP(v : \)=J2bj(v,l)\ j . 

3=0 

In |Gal| , the first author showed that the colored Jones function Jk of a knot K satisfies an essentially 
unique smallest degree g-diffcrcncc equation PkJk = where the coefficients a,j(u,v) of Pk are rational 
functions of u and v with rational coefficients. 

In |Gal| . the operator Pk was called the non- commutative A-polynomial of K . 

In |Gal| . it was conjectured that: 

Conjecture 1. (A3 Conjecture) Up to a multiplication by a polynomial in v, we have 

chP K (X,v) = A K (L, M)\( Lt M*)=(\,v) 
where Ak is the A-polynomial of K , defined by CC GLS) . 

The A-polynomial of K parametrizes the moduli space of characters of SL2(C) representations of tti(S z — 
K), restricted to the boundary torus dM. The A-polynomial of a knot is an important ingredient to the 
Geometrization of the knot complement and its Dchn fillings. 

The A-polynomial Ak of a knot K in S 3 satisfies symmetries, which we will list here, and refer to [HCGLS 
and |GLj for proofs. 

(51) It has integer coefficients and even powers of M, that is Ak{L, M) E Z[L, M 2 ]. 

(52) It is reciprocal, that is, A K {L- l ,M- x ) = ±L k M l A(L,M). 

(53) It is tempered, that is the edge polynomials of its Newton polygon are cyclotomic. 

(54) It specializes to 

A K (L,l) = ±(L-ir+(L + l) n - 

for some integers n±. 

(55) L — 1 is always a factor of Ak, that corresponds to U(l) representations. 

If the colored Jones function of a knot was an iS-regular solution to a regular g-diffcrence equation, and 
if the AJ Conjecture were true, then it follows that for every a £ [0, 1] we have: 

(45) lim '1 = ai K {a) = ai K (a) 
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where <jg K is the A-entropy of a knot, defined as follows. 

Definition 7.1. For a knot K in S 3 , let Lj(t) for j = 1, . . . , d denote the roots of the equation 

A K (L j (t),e it / 2 ) = 

for t € [0,27r], where d is the L-degree of Ak- Fix a partition U pe -p/ p of [0, 2ix] by closed intervals with 
nonoverlapping interiors and a permutation o~ p of the set {1, . . . , d} such that 

\L ap{1) (t)\>\L ap{2) (t)\>--->\L ap{d) (t)\ for all t e I p . 

For every collection S — {S p \p G V} of subsets of {1, . . . , d}, we can define the S- entropy 

a£ K : [0, 1] — > R 

by 

1 f 2 ^ 

^S,k{°) = TT logxs(at)dt, 
27r Jo 

where \s ■ [0, 1] — > R is defined by 

Xs(t) = max|L CTp(i) (i)| if t6J p . 

It is natural to ask how the A-entropy of a hyperbolic knot (evaluated at a = 1) compares to the 
Hyperbolic Volume. The answer to this question is essentially contained in work of D. Boyd, |Bo| . which we 
quote without proof here. We urge the reader to look in |Bo| for beautiful and suggestive calculations. 

Boyd introduced and studied another invariant of a knot, the Mahler measure 

m K = / log\A K (x,y)\dxdy 

JS 1 xS 1 

Using Jensen's formula, and the symmetry (S3), it follows that 

d „2tt 

3=0 Jo 

where log + (a) = logmax{l,a}. 

Using (S2), it follows that \/Lj(t) is an eigenvalue for every eigenvalue Lj(t). 

More generally, among the roots Lj(t) there is a distinguished one, L\, corresponding to the discrete 
faithful representation when t — 2ir. Let Ld{t) — 1 denote the eigenvalue corresponding to the U(l) 
representations. Boyd informs us that for 2-bridge knots K (in particular, for the 3i and 4i knots), it is true 
that 

a {l,d}.K C 1 ) = Y0l K- 

It follows by (S4) that the eigenvalues collide at t = 0. Moreover, Boyd informs us that for 2-bridge knots 
there exists a to G (0, ir) such that |£j(i)| = 1 for all j and all t : to < t < it. 

Thus, if we want to apply Theoremn]to the GHVC, we need to deal with irregular ^-difference equations. 
We will discuss this topic in detail in a later publication. 

Meanwhile, let us discuss some examples, taken from |Gal) . 

7.2. Examples: The 3i and 4i knots. In this section we discuss in detail g-difference equation of the the 
colored Jones function of the two simplest knots, namely the trefoil 3i and the figure eight 4 X . The former 
is not hyperbolic, and the latter is. 

In |Gal| . the first author computed that the colored Jones function (resp. J4J satisfies the second 
(resp. third) order g-difference equation 

Pi! 'hi = resp. P 4l J 4l = 
where the noncommutative A-polynomials Ps 1 and P 4l are given by: 
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ft, = 



g 3 Q 2 (g 2 ~g 2 Q) 

q 3 - q 4 Q 2 

(q - q 2 Q)(q + q 2 Q){q A - q 5 Q + q 6 Q 2 - q 7 Q 2 - q 7 Q 3 + q a Q 4 ) 



q 2 Q(q~q 4 Q 2 )(q 3 -q 4 Q 2 



E 



-l + q 2 Q _ E 2 



4i 



Q(q~q 4 Q 2 )' 
q 5 Q(-q 3 + q 3 Q) 
(q 2 + q 3 Q)(-q 5 + q 6 Q 2 ) 
(q 2 - q 3 Q)(q 8 - 2g 9 Q + q 1Q Q - q 9 Q 2 + q W Q 2 ~ q n Q 2 + g 10 Q 3 - 2g n Q 3 + q 12 Q 4 ) 

q 5 Q(q + q 3 Q)(q 5 ~q 6 Q 2 ) 

'0 2 -^n w 3 4-n 11 3 + n 12 4 ] 



E 



(-<? + q 3 Q)(q 4 + q 5 Q - 2<z 6 Q - q 7 Q 2 + q s Q 2 - q 9 Q 2 - 2 q w Q 3 + q n Q 3 + q 12 Q 4 ) 2 



q 4 Q{q 2 + q 3 Q)(-q + q & Q 2 ) 



q 4 Q(-l + q 3 Q) e3 



(q + q 3 Q)(q - q e Q 2 ) 

If we wish, we may clear denominators in P 3l and P 4l . It follows that the characteristic polynomials are 
given by: 



chP 4l {L,M) 



M(l + M) 

(i - 1){L - LM - M 2 - 2LM 2 - L 2 M 2 - LAI 3 + LM 4 ) 



M(l + M) 2 

Inspection shows that P 3l and P 4l are not regular. Nevertheless, let us try to compute the S'-entropy. 
For the case of 3i, we have \Lj{t)\ = 1 for j = 1, 2, 3 and in this case 

4 3l (l)=V l 3l =0 

for all S. 

For 4i knot, we have 3 eigenvalues L\{t), L2(t) = l/Li(t) and L 3 (t) — 1. Assuming appropriate choices 
for the branches of the eigenvalues, the plot of log |Zi(t)| and log \L2(t)\ = — log (£) | fori G [0,2tt] is given 
by: 




It follows that 

«rf 1)3})4l (l) = vol(4 x ) = 2.029883 ofi,^},^ 1 ) = 2vol 4i = 4.05977. 

Since the HVC is true for the 4i knot, it suggests that the colored Jones function lies in a strictly smaller 
subspace of the vector space of solutions to the ^-difference equation PrJk = 0. Using work of Murakami 
|Mu| . one can figure out exactly the selection principle; that is which locally fundamental solutions contribute 
to the colored Jones function. 
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Note that the associated g-difference equation of the 4i knot has the following features: of collision, 
resonance and vanishing: 

• The eigenvalues collide at t — (since Li(0) = Li(0) = —1), at t = ir/3 (since Li(n/3) — Li(tt/3) = 
— 1) and by symmetry at t — 27r/3 and t = 2tt. 

• There is resonance on the interval [tt/3, 2ir/3] where all three eigenvalues have equal magnitude. 

• There is vanishing of the coefficients at t = ir/2 (since the denominator M + 1 of the coefficients is 
singular at t — n/2). 

Moreover, there is an additional difficult problem of selection principle. 
We plan to study these problems in later publications. 
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